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Abstract 

We show that the Bronnikov theorem on the nonexistence of regular 
electrically charged black holes can be circumvented. In the frame 
of nonlinear electrodynamics, we present the exact regular black hole 
solutions of a hybrid type. They are electrically charged, but contain 
a 'dual' core confining a polarization of magnetic charges. 

The considered example is based on a modification of the Ayon- 
Beato &; Garcia solution. It represents a very specific realization of 
the idea on confinement based on dual electrodynamics and dual su- 
perconductivity. 
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1 Introduction 



Last years considerable interest in regular black hole solutions has been re- 
newed |], H, HQ, and in particular, to those based on nonlinear electrody- 
namics (NED) j|, |5|, |6], [7], |8j . Analyzing the Ayon-Beato & Garcia black hole 
solutions from nonlinear electrodynamics j|, |5|, [|, Bronnikov presented three 
important results 0: 

(i) - a version of no-go theorem claiming the nonexistence of the regular 
electrically charged black holes and particlelike solitons with regular center 

ft 

(ii) - the existence of cusps and branches in the Lagrangian used for the 
Ayon-Beato - Garcia solution ||; and 

(iii) - a class of exact regular magnetically charged black holes and soli- 
tonic solutions. 

He showed that the branch, realized in the Ayon-Beato - Garcia solution 
near the core, does not tend in the weak field limit to the Maxwell theory 
and presented two theorems leading to the conclusion that regular black hole 
and solitonic solutions can be only magnetically charged. 

The unusual properties of the Ayon-Beato - Garcia solution obtained by 
Bronnikov have captured our attention in connection with recent suggestions 
on a phase transition which has to occur for regular black hole and parti- 
clelike solutions in the core region |], ||, |3|. In this work we show that the 
conditions of the Bronnikov theorem can be circumvented by dealing with 
the models having a phase transition near the core of solution, when electric 
field does not extend to the center of the solution. 

After some modification of the Ayon-Beato - Garcia solutions, we ob- 
tain a new type of regular, electrically charged black hole solutions. The 
main peculiarity of these solutions is a 'magnetic' core described by dual 
electrodynamics. This class of solutions deserves interest not only as an ex- 
plicit example of the regular electrically charged black hole solution, but 
also as an example of particlelike solution with a very specific realization of 
the ideas on quark confinement based on dual electrodynamics M and dual 



superconductivity [QG, 11, 12, 13, 14 



1 See also list of references in jjj. 
2 As well as wormholes or horns. 
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2 NED and F-P Dual representations 

The action for NED in general relativity is 

S=^Jd 4 x^[R-C(F)}, (1) 

where R is the scalar curvature, F = F tlv F tlv , F^ u = d^A u — d u A^, and £(F) 
is a nonlinear function. 

The equations, following from Eq. ([[]) are 

V fl (L F F^) = 0, V M *i^ = 0. (2) 

where L F = dC/dF, and the equations V M * F^ u = are the Bianchi identi- 
ties. 

By introducing the electric intensity E and magnetic induction B, 

E = {F l0 }, B = = 1^°%,, (3) 

one can obtain from Bianchi identities V^-kF^ = the pair of Maxwell equa- 
tions VB = 0; V x E = -dB/dt. The dynamic equations V^LpF^) = 
can be expressed via the tensor 

P» v = dC/dF^ = L F F^, (4) 

in the form 

V M P^ = 0; (5) 
that, in the terms of the electric induction D and magnetic intensity H, 

D = {P l0 } = L F {F l0 }, H = {*P° 1 } = L F {-kF 0i } : (6) 

yields the second pair of the Maxwell equations VD = 0; V x H = 
dD/dt. These equations show that in NED L F plays the role of electric 
susceptibility e = L F and the magnetic permeability is given by \i = l/L F . 

Relation (HI) allows one to rewrite the field equations (0) in an equivalent 
dual P form 

V,(?n = 0, V^*P^ = 0. (7) 

L>F 

3 We use signature — V ++■ The Hodge star operator is given by *F A " y = \rf va ^ F a $, 
Vftvap is completely skew-symmetric, 7/0123 = 
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It should be mentioned that it is a dual description of the same physical sys- 
tem. This description can be obtained from the electromagnetic Lagrangian 
T~t(P) ( expressed via P = P^ U P^ U = L 2 F F ) determined by Legendre trans- 
formation 0, 0, [L5|] 

H{P) = 2FL F - C. (8) 
The following relations can be easily obtained 

C = 2PH P -H] H P = dH/dP = 1/L F . (9) 

The Hodge star duality operation allows one to give a description of the field 
intensities via dual components F^ = -kF^ = ~r] "PF a p. Since *★ = —1, by 
using the Hodge dual description the following relations hold: 

F = -F- P = -P; (10) 

and also 

L F = -L F ; H P = -Hp- L(F) = L(-F); H(P)=H(-P); 

(11) 

and the equations @ take in this description the form 

V M (L F * F"") = 0, V M F^ = 0. (12) 

Note that the Hodge dual description changes sign of Lp that can lead to 
nonphysical values of electric susceptibility e < 0. 

Vice versa, if the Lagrangian leads to Lp < it can be considered as a 
sign to pass to another, physically valid (Hodge) dual description. 

3 Nonexistence of Electrically Charged Reg- 
ular Black Holes? 

The Bronnikov no-go theorem is based on the assumption that the NED 
Lagrangian has to lead to Maxwell theory at small F: C(F) — > F as F —>■ 
0. This is equal to Lp — > 1 at small F. He shows that real behavior of 
electrically charged solution near the center leads to finite values of \FLp\ 

4 See, for example, [|l6] 
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and to divergence of \F\L 2 F . Therefore at the center the ratio FL ^ = Lp 
tends to infinity that is non-Maxwellian behavior. 

This derivation is correct. However, there are two points which throw 
doubts on the initial postulates. The first one is connected with the assump- 
tion that the electromagnetic field is extended up to the center of the solution. 
At least in flat space there is a remarkable example of the electrically charged 
solution, present only in the region r > a. It is the classical Dirac electron 
model, as a charged sphere. (] The second argument is connected with the 
demands of Maxwellian behavior. Indeed, the condition 1 as F ^ 

(to recover Maxwell theory) is reasonable at spatial infinity. However, at 
the core, although the gravitational field is regular everywhere and becomes 
weaker as one approaches the center (the metric tends to flat space-time), 
the spatial volume is much smaller (on many orders) than in the asymptotic 
limit where gravity is again weak and where Lp — > 1 makes sense. Thus 
both types of boundaries are essentially different. Another argument can be 
a strong vacuum polarization and the presence of other ( maybe virtual quan- 
tum ) fields which have a straight relation to the formation the nonvacuum 
values of Lp = e. Therefore C(F) can have a nonlocal behavior depending 
not only on the value of F but also on full field configuration and proximity 
to the core. 

If one restricts oneself to Lagrangian densities which are local functions 
of F, then, Bronnikov is completely right. But assigning essentially the same 
physics to the core region as to the asymptotic region (in electromagnetic 
terms) is not convincing enough. Therefore his restrictions for regular so- 
lutions are too especial to be accepted in general. This point of view is 
suggested in fact by the given Bronnikov analysis showing the existence of 
cusps and branches in the Lagrangian of the Ayon-Beato - Garcia solution. 
In fact, his analysis shows that the Ayon-Beato - Garcia Lagrangian has 
different behavior at the core and far outside the body, although F — > in 
both cases. 



5 The Bronnikov theorem 2 tries to exclude this case showing nonexistence of a horn 
in space-time. However, there can be other (nonmetrical) mechanisms restricting the 
penetration of the electromagnetic field into region r < a, in particular, superconducting 
properties of core or a special behavior in NED as it will be shown further. 
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4 The Ayon-Beato — Garcia Solution 



The electrically charged regular black hole solution has been given by Ayon- 
Beato - Garcia in [||]. Most investigations in nonlinear electrodynamics are 
connected with the Born-Infeld type of nonlinearity (see, e.g., |L7|, [18], |19| 
and references therein) . However, as it was mentioned for regular black hole 
solutions in 0, the "...Born-Infeld Lagrangian is useless in this topic, since 
it gives rise to a singular black hole solution, at least in the case of spherical 
symmetry pi|..." , and Ayon-Beato - Garcia start from the nonlinear function 



H(P) = P/cosh 2 {s(-2g 2 P) 1 / 4 } 



(13) 



where s = \q\/2m (q is the electric charge and m is the mass of the body). 




Figure 1: Dependence C(F) for Ayon-Beato - Garcia solution. 

Note, that particular form of nonlinearity 7~t(P) is not essential, and 
existence of regular black hole solutions is determined by the asymptotic 
properties of this function. What is important is the Maxwellian behavior 
on large distances, 7i(P) = P at small P, and a rapid fall-off of this function 
by P — > oo. These demands lead to the existence of an extremum of the 
function TC, and moreover to the existence of at least one sector where Hp 
is negative. For instance, we use for illustration the simplest function of this 
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sort 

H(P) = Pexp{-\P\}. (14) 

The corresponding Lagrangian in L description is C(F) = 2HpP — 7i(P), 
where P has to be expressed via F = P / L F , ( see Fig. 1 and upper graph in 
Fig. 2.) . 

It yields the system of electromagnetic equations 

V,(Pn = 0; (15) 

and Einstein equations = where the energy-momentum tensor is given 
by 

T;= l -{AL F F, a F^~5;C). (16) 

In the case of spherical symmetry the electrically charged solution of Eq. 
©is 

Por = -Po = q/r 2 (17) 

or 

F 0r = -F r0 = L~ F l q/r 2 - (18) 

the other components being equal to zero. Thus P = —2q/r 4 . The Einstein 
equations can easily be solved leading to the mass term M(r) = m[l — 
tanh(g 2 /2mr)] and to the metric 

ds 2 = -[1 - 2M(r)/r]dt 2 + [1 - 2M(r)/r]- l dr 2 + r 2 dtl 2 , (19) 

which is regular by r going to zero. 

As mentioned by Bronnikov [^], the Lagrangian given by the smooth and 
regular function TC(P) displays cusps and branches in F description. Appear- 
ance of cusps at the points Pi and P3 is connected with multivaluedness of 
the inverse map P(F) ( see fig. 2) and is a general feature of the considered 
class of nonlinear functions. Very specific behavior of function L(F) by large 
P, on the section P 2 — P3 — P4 is determined by the following phenomenon: 
The magnetic permeability /i = Hp is going to zero as P — > P2, the point 
of TC extremum. Consequently, electric susceptibility e = 1/Hp is divergent 
at this point. The sphere, determined by the equation P = P 2 looks like an 
ideal conducting surface, and electromagnetic field F^ ( and F) will tend 
to zero approaching this surface. The solution, analytically extended inside 
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Figure 2: Graphs of TL(P) and F(P) for Ayon-Beato - Garcia solution and mod- 
ified solution F mo d(P). 



this sphere departs drastically from a Maxwellian behavior. Nevertheless it 
still has reasonable properties in many respects. For instance, the electro- 
magnetic field F^ u and the invariant F go both to zero as the very center of 
the solution is approached ( although L F — > oo at P 2 and P 4 ) and the metric 
is regular everywhere, leading to regular black hole solutions . 

In spite of these properties, the physics inside the core looks rather un- 
usual since the value of electric susceptibility e = 1/ 'Hp is negative there and 



does not correspond to any known physical mediums [pl|. Combining NED 
and general relativity, one can show from the system of equations (|l5l), (|16|) 
that the corresponding contribution of effective charges to the metric will be 
of opposite sign in the regions with negative Lp, as if the charge would be 
imaginary valued there. 



5 Modification of the Ayon-Beato — Garcia 
Solution 

The fact that regularity of black hole solutions can be reached without 
troubles in the magnetically charged solutions suggests the idea to combine 
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the electrically and magnetically charged solutions by using a phase transition 
near the core to magnetic phase, regularizing the black hole singularity. In 
this case, the existence of different branches in the Ayon-Beato - Garcia 
solution represents a perfect possibility for realization of this idea by means 
of matching external and internal solutions with necessary properties. As 
we shall see, the transfer to dual electromagnetic phase near the core of 
the Ayon-Beato - Garcia solution allows also one to avoid abnormal branch 
(P 2 — P3 — Pa) with negative L F . 

Looking on the graph L(F) (Fig. 1) one sees that it can only be achieved 
by replacing this branch by the branch obtained by a mirror map with respect 
to the L axis, as it is shown in Fig. 3. This replacement takes place in the 
region of strong field \P\ > IP2I, i-e., inside of the sphere r < r = Ifr] 1 / 4 
surrounding the singular point of P(r). 



-L 




Figure 3: Initiate and modified dependencies L(F). 



As we shall see this replacement corresponds to "dual" electrodynamics 
with a polarization of magnetic charges inside the sphere r = r , while the 
solution for r > r survives to be the electrically charged Ayon-Beato - 
Garcia one. Matching of the modified internal solution and the external 
Ayon-Beato - Garcia one is smooth with respect to the metric and physical 
field intensities. 
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The properties of modified Lagrangian L(F) and the modified relations 
will be the following. 

1. The region of definition of Lagrangian L is extended to positive values 
of the invariant F = B 2 — E 2 , and modified electric susceptibility, 
e = Lp = —Lp, is positive. 

2. Matching of the new branch with the initial one is smooth and occurs 
at the point P 2 where F = F = and L(F)\p 2 = L(F)\p 2 . 

3. The modified dual intensities are related as 

P^ = LpF^; P = {LpfF. (20) 

At the point P 2 electric susceptibility is divergent Lp\p 2 = L F \p 2 = 
+oo, while the magnetic permeability /i = 1/Lp is bounded and C . 

4. As far as in the modified region F > we get also P > 0. On the other 
hand P = H 2 — D 2 is negative in the region r > tq (as in the case of 
an electrically charged solution). Therefore the modified function P(r) 
will be discontinuous. The change of sign of the invariant P at the 
point P 2 is evidence of a transfer to magnetic phase by r < r . 

Due to the Hodge star operation identity ** = — 1, the equality P = —P 
can be resolved as P^ u = ± * P^ v . Whence, the relation (^j shows that the 
modified Ayon-Beato - Garcia solution (|l"8D by r < r describes a field of 
magnetic induction in radial direction n 

B = TqnLf/r 2 . (21) 

Magnetic induction B characterizes the true intensity of the field strength 
in media (T^|. In modified solution it tends to zero at r = and at r = tq. 
Contrary to the standard Maxwell equation VB = describing the absence 
of magnetic charges, we obtain here the magnetic charge distribution 

dF 

VB = J mag = T(q/r 2 )nVL F = T ( q /r 2 )L FF —, (22) 

displaying realization of dual electrodynamics in the core region and singular 
concentration of magnetic charges at r = and r = r$. 

The initial and modified dependencies TC(P), F(P), and L(F) are given 
in figs. 1 and 3. 
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Figure 4: Radial dependencies of L and F for modified solution. 



5.1 Stress-energy tensor and metric. 

We use the Kerr-Schild form of metric [|T], [] In the case of spherical 
symmetry 

9^u = V^ + [f(r)/r 2 }k^ (23) 

where rj^ is the auxiliary Minkowski metric and = {l,n}. Co- and 
contravariant forms of the metric in the Kerr angular coordinates {t, r, 8, 0} 
are given in the Appendix. We have 

F 01 = -F 01 ; F 23 = r 2 sm6F 01 ; F 23 = F 23 / (r 4 sin 2 9); (24) 

which, taking into account Eq. (|T8|) , yield 

F = 2F 23 F 23 = 2(F 01 ) 2 = 2L F 2 q 2 /r\ (25) 

Similarly, 

P = 2P 23 P 23 = 2(P i) 2 = 2q 2 /r 4 . (26) 

6 See also where it was shown that the static solutions of the Einstein equations 
with spherical symmetry and Tq = T± belong to the Kerr-Schild class. 

In this section we use slightly different units, adapting to the works |2^, [| and setting 
G = 1. The reason for this is the possibility of using the Kerr-Schild formalism for further 
extension to rotating solutions. 
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The modified energy-momentum tensor inside the core acquires the form 
T M„ = J_^~ L _ A p^p^/~ Lp) = -±diag{L, L,L- 2P/L F , L - 2P/L F }. 

Outside the core, by r > r it is 



(27) 



T^> = -±(5»L-4:P m P™/L F ) = -±dzag{L-2P/L F ,L-2P/L F ,L,L}, 

(28) 

where P = —P = —2q 2 /r 4 . As far as contributions of terms 2P/Lp and 
2P/L F , together with their first derivative, tend to zero approaching r , one 
sees that these expressions match smoothly at r . 



- L , H 




0.5 



Figure 5: Radial dependencies of — L and H. 

Comparing corresponding expressions for the generalized Kerr-Schild class 
of metrics Jl] 

T; = {9m)- x diag{-2G, -2G, D + 2G, D + 2G}, (29) 
where (see the Appendix) 

D = -f"/r 2 , G=(f/r)'/r 2 (30) 
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and ' = d/dr gives the relations: 
For external field 



if fr)' 



Mr), 



(31) 



2vrr 2 



where p(r) is the energy density in this region. 
For the external part of the function f(r) 



f(r) = rM (ext) (r), M (ext) (y) 



2vr [ V H (ext) (r)r 2 dr + C, (32) 



where C is the integration constant. 

For function f(r) inside the core, r < r , we obtain 




(33) 



that gives rise for the internal region to 



f(r) = rM {int) (r), M (int) (y) = 2ir F L(r)r 2 dr. 

<J 



(34) 



Each of the functions M^ ext / int ){y) plays the role of the mass confined inside 
the sphere y < r. Constant C can be determined by matching both sectors 
of f(r) at r = r that yields C = M( int )(r ). This matching is smooth as far 
as TC(r ) = —L(r ), thus the metric (|23|) is smooth everywhere. 

By analyzing dependencies L(r) and Ti.(r), Figs. 4 and 5., one sees that 
there is a core region (r < 0.6ro) which does not give contribution to mass 
at all. The space is flat in this core since f(r) ~ there. The region 
0.6 < r/Yo < 1 is a transitional one with a dual magnetic phase, whereas the 
region 1 < r/r < 1.6 is one with an electric phase. In the region r > 1.6ro 
we have Lp ~ 1, it is practically a Coulomb phase. 

Let us estimate characteristic scale r for the Ayon-Beato - Garcia type 
of nonlinearity flTBD . Radius r is determined by the equation Hp = 0, or by 
the root of equation £tanh(£) = 2. It yields ( = 2.066, and 




(35) 



where 



r e = q 2 /2m is a classical electromagnetic radius. 
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For comparison, in the Born-Infeld theory (in our notations), Lagrangian 

is 

Cbi = ^{VT+WF -1), (36) 
and has a dual Lagrangian 



Hbi = ^(1 - Vl-26 2 P). (37) 

The typical NED branches are absent here, and the sphere Hp = shrinks to 
the point r = 0. Therefore the considered above type of solutions is absent 
in Born-Infeld theory. 

Black holes in Born-Infeld theory were actively studied (see |T3fl which 



includes the original references). The Born-Infeld black hole solution has 
finite total energy E = m w 1.236e 2 /ro and a smooth distribution of effective 



charge p e ff = -^div'E in the region ~ r = yeb . Parameter b~ l = e/rg 
characterizes a critical value of the electromagnetic field. 

Considering particlelike solutions, one sees that, similarly to NED, the 
effective radius of the elementary particle in Born-Infeld theory is close to 
classical electromagnetic one 

r = 2.472r e ~ i(T 13 cm. (38) 

However, in contrast to NED, the Born-Infeld solution has a 5-function source 
at the origin, leading to a conical singularity in curvature. Therefore the 
Born-Infeld theory, as well as NED, has classical radius r e as one of the 
characteristic scales. However, regularizing ( or smearing ) the black hole 
singularities, both the theories display that their minimal ( space ) scale is 
extended up to Planck sizes. 

On the other hand, the Born-Infeld theory can be derived from open 
string theory |TJ| |TB|, and parameter b takes in this case the value b = 
2 r na' = T -1 defined by Regge slope a'. On the level of particlelike solutions, 
it does not contradict with b^ 1 ~ e/rg, leading to approximate value of the 
fine structure constant a ~ ( see, e.g, Hl8|). However, in the modern 
superstring theory tension T( ss ) defines the critical scale for all field strength, 
including curvature, and characteristic radius r crit = Joi'^ corresponds to 
the Planck scale. It shows that the scale range of nonlinearity in NED action 
is much larger than the Planck scale that governs higher order corrections in 
curvature. 



14 



6 Concluding remarks 



The Ayon-Beato - Garcia solution contains distribution of electric charges 
JY e s = V^F^ near the sphere r = r$ and in the center. In the presented modi- 
fied solution the charges inside the sphere r = r are magnetic JY m s = V ^F^ 
fulfilling the Dirac idea on dual electrodynamics 0. These charges are con- 
fined inside the sphere r = ro, forming a magnetic vacuum polarization. For 
the external observer this solution exhibits only electric charge. Since vac- 
uum of the core expels electric charges one can speculate that it possesses 
superconducting properties, whereas the vacuum of the external observer 
demonstrates the dual superconductivity, expelling the magnetic charges. 
The considered matching of the electric and magnetic phases corresponds to 
a sharp phase transition between these vacua ( London limit of superconduc- 
tivity). 

The similar baglike structures with a phase transition to dual supercon- 
ductivity ( approximate solutions ) were recently discussed for regular rotat- 
ing black holes to supergravity 0. It was assumed that superconductivity 
is caused by the Higgs ( chiral ) fields forming two dual ( supersymmetric ) 
vacuum states and a smooth domain wall interpolating between them. It is 
remarkable that nonlinear electrodynamics allows one to get exact and self- 
consistent solutions of this sort belonging to the Kerr-Schild class of metrics. 
However, it could seem wonderful since the Higgs fields providing the vacuum 
states and phase transition are absent here. On the other hand, the NED 
Lagrangian L(F) can be represented in the form (^p-)F showing that the 
factor ^p- can be identified as a dilaton ( corresponding to a metric on the 
moduli space in supergravity). It allows one to suppose that NED can be 
assumed as a limiting case of supergravity when the potential is very sharp 
and chiral fields are dropped out. The presented exact solution from NED 
shows that dilaton has to be very essential in the baglike models. 

A parallelism between the above baglike models and the field models of 
dual superconductivity should also be mentioned. Indeed, there are known 
the dual superconductivity models of two type: the models with Higgs fields, 
leading to solutions with a soft phase transition ( finite thickness of tran- 
sitional layer ), see, e.g., [12|, [L3|], and the similar to NED models without 



Higgs fields, see, e.g., |TT|, [L4|], demonstrating the sharp phase transition con- 



nected with the London limit of superconductivity. The field models of the 
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first type represent in fact a truncation of the four-dimensional supergravity, 
while the models of the second type represent, apparently, its singular limit. 
[] However, these questions demand more detailed consideration elsewhere. 

It would also be very interesting to look on NED from a high dimen- 
sional point of view, getting it as a result of compactification or as a four 
-dimensional slice of a multidimensional spice-time in the spirit of early brane 
world work 



2J. 



Finally, taking into account nontrivial topology of the Kerr geometry, its 
exceptional character with respect to dual rotations [15 and relationships to 
the spinning particlelike solutions ( see, e.g., [55] and references therein ), 
the treatment of corresponding regular rotating solutions is very important 
and could bring some surprising results. In particular, the appearance of a 
closed vortex string on the place of the Kerr singular ring is expected. 
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Appendix. The Kerr-Schild form of metric in 
the case of spherical symmetry 

The Kerr-Schild form of metric [22, II] in the case of spherical symmetry is 

9i*v = Vijv + Ihk^kv, (39) 

where rj^ is the auxiliary Minkowski metric, h = f{r)/r 2 , and k^ = {l,n}. 
In the Kerr angular coordinates {t, r, 9, 4>} we have 







„2 • 2, 



/ 2h - 1 2h 

2h 1 + 2h 



V r 2 sin^ J 



(40) 



7 Note that effective Higgs fields can also be introduced in the models of the second 
type p). 
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where h(r) = f{r)/r 2 , \f—g = r 2 sin#, and the contravariant form of metric 
is 



( -(1 + 2/i) 2h 
2h l-2h 
1/r 2 









V 



(r 2 sin 2 i 



1 I 



(41) 



The energy-momentum tensor for the generalized Kerr-Schild class of metrics 
is 1,| 

T; = (87r)~ 1 dm(/{-2G', -2G, D + 2G, D + 2G}, (42) 

where 

D = -f"/r 2 , G=(f/r)'/r 2 , (43) 
the prima means d/dr, and energy density is p = 2G/87T. 
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